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Understanding the nature of low energy quasiparticle 
(QP) excitations in isolated vortices and vortex lattices 
of type-II superconductors is of considerable importance 
since for low temperatures T these excitations impact on 
various static and dynamic properties. The flux lines are 
known to affect the low energy QP states through the 
following mechanisms: (i) QP scattering on the gap in- 
homogcnuity in vortex cores (which is known to result in 
formation of localized core states in s-wave compounds) 
; (ii) QP scattering on the potential proportional to the 
superfluid velocity Vs (iii) long-range magnetic field 
effects which are responsible for a finite curvature of a 
quasiclassical trajectory in a vortex lattice (iv) inter- 
action with a pure gauge potential due to the Aharonov- 
Bohm (AB) scenario [@,||. For conventional s-wave su- 
perconductors these mechanisms have been studied for 
several decades and now one may conclude that the phys- 
ical picture of the electronic structure of the mixed state 
in s-wave systems is rather clear. The situation is dra- 
matically different in superconductors with gap nodes 
where the vanishing pair potential in the nodal directions 
results in qualitative changes in quantum mechanical mo- 
tion of QPs. The interest to these fundamental issues is 
stimulated by recent experimental observations of uncon- 
ventional behavior of QP excitations in the mixed state 
of high-Tc compounds [|5|"pT|] where the dominating order 
parameter is believed to be of d-wave symmetry. 

In this Letter we focus on the quantum mechanical ef- 
fects caused by the interaction of QPs with pure gauge 
potentials induced by flux lines. The main goal of our 
theoretical analysis is to show that contrary to the s- 
wave case the AB effect plays a crucial role in the be- 
havior of low- lying QP states in d-wave compounds. Let 
us consider a single isolated vortex line which carries the 
flux quantum 0o — ■KTic/\e\. Due to the Meissner ef- 
fect a magnetic field H is screened at the London pen- 
etration depth Al and at large distances from the vor- 
tex center r 3> the vector potential takes the form: 
A = (/)o[zo, r]/(27rr^), where [zo,r] is a vector product of 
vectors zq and r, and Zq is a unit vector chosen along 
the vortex axis. Such a vector potential can not be ex- 
cluded from the Bogolubov - de Gennes (BdG) equations 
(describing the quantum mechanics of QPs) using any 
single- valued gauge transformation. Note that within the 
standard Schroedinger theory for a particle with a charge 
e an AB solenoid with a magnetic flux $ = (/jq is known 
to be the most effective scatterer [|l2|. However one can 
see that for s-wave superconductors the effect of the AB 
potentials on the QP density of states (DOS) should be 
rather small. The point is that the size of the semiclas- 
sical wavepacket propagating outside the flux tube (i.e., 
in the region r 3> \l) appears to be much less than 
the impact parameter. As a result, the interference ef- 
fects are small and can be neglected. As we see below, 
such a conclusion is no more valid if we consider super- 
conductors with gap nodes. Hereafter we assume Fermi 



surface (FS) to be two-dimensional (2D), which is ap- 
propriate to high- Tc superconductors, and take the gap 
function in the form = 2Aokxky/kp (the x axis is 
taken along the [110] crystal direction and thus makes 
an angle 7r/4 with the a axis of the Cu02 planes). The 
nature of low energy QP excitations in the vortex state of 
d-wave systems attracted recently a great deal of atten- 
tion (see pP,p^-p6[| and references therein). In contrast 
to conventional superconductors, the DOS at low energies 
s in d-wave systems is dominated by contributions which 
come from the regions far from the vortex cores and 
associated with extended QP states with momenta close 
to the nodal directions. This conclusion based on the 
semiclassical approach has been confirmed by the recent 
numerical analysis of the BdG equations for a single 
isolated vortex line (in the limit A^ — + oo). Note that 
the calculations presented in [|l9| also point to the ab- 
sence of truly localized core states or any resonant levels 
in the pure d-wave case though such states were observed 
in numerical simulations |l8| . 

Let us start with some qualitative arguments which 
indicate a significance of pure gauge potentials in the 
quantum mechanics of extended QP excitations. In the 
homogeneous state the low energy excitations (confined 
to one of the gap nodes) are described by a Dirac-like 

spectrum: = ±fi^Vpq\ + ' where Vp is the 

Fermi velocity, and {q\\,qi_) defines a coordinate system 
whose origin is at the node, with q±^ {q\\ ) normal (tangen- 
tial) to the FS. The Va value determines the gap slope 
near the node and for our choice of the gap function is 
given by the expression Va = 2Ao/(?ifcF). One can sep- 
arate the following length scales for QP wavefunctions: 
an atomic length scale kp^ and two characteristic wave- 
lengths of a slowly varying envelope Zj^ ~ gj^ ~ KVp/e, 
h ^ HVa/s- The length scales l± and /y deter- 

mine the size of the semiclassical wave packet, which 
appears to diverge in the low energy limit. Such a di- 
vergence is responsible for a crucial role of quantum me- 
chanical effects in QP motion and, in particular, for the 
extremely large AB scattering amplitude. This conclu- 
sion is in a good agreement with the recent calculations of 
the QP band spectrum P, p2| , p4[[25| ] and thermal conduc- 
tivity |2^ in the regime Hd <C <C Hc2- Our further 
consideration of the low-energy QP states (|e| <C Aq) 
with momenta locked to the nodes is based on the essen- 
tially 2D quantum mechanical model proposed by Simon 
and Lee [0. Let us introduce the gap function describ- 
ing the mixed state of a d-wave superconductor in the 
form A(k, R) = Ad(k)vI/(R), where ^'(R) = fexp{ix) is 
the order parameter used in the Ginzburg-Landau the- 
ory. In the limit fcj?^ ^ 1 (^ is the coherence length) 
one can divide out the fast oscillations on a scale kp^ 
in the particle-like and hole-like parts of the wavefunc- 
tion (u^v) = {u,v)exp{i'kFv) and simplify the nonlocal 
off-diagonal terms in BdG equations. Omitting details 
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of the derivation, we start from the BdG equations hn- 
earized in gradient terms for quasiparticles with momenta 
close to a certain gap node direction (e.g., ki — {kp, 0)): 
Hg = eg, where g = {uexp(—ix),v), and the Hamilto- 
nian H outside the vortex cores takes the form: 



H = Ho + \e\ip 



(1) 



9 Xy^ 



Here ax,^z are the Pauli matrices, the scalar poten- 
tial (p = MVfVsx/\g\ results in the Doppler shift of 
the QP energy, M is the electron effective mass, Vsx 
and Vsy are the components of the superfluid velocity 
Vs = FsxXo -I- Vsy-ya, xq, yo, zq are the unit vectors 
of the coordinate system with Zq chosen along the c- 
axis. For an isotropic FS such a description is correct 
only for e <C Ao/(fc;-^), but the range of validity may 
become larger if the FS is somewhat flattened at the 
nodes. The single-valued transformation of wavefunc- 
tions {g = {uexp{—ix), v)) used above was introduced 
previously in [ p2|j23| . Note that other possible versions 
of single-valued gauge transformations were discussed in 
1^,^. The superfluid velocity Vs(r) can be written as 
a superposition of contributions from individual vortices 
situated at points r^: 



V.(r) 



h \ ^ [zo,r- 
^ Ailr- 



2M 



where Ki is the Mcdonald function (modified Bessel func- 
tion) of the first order. Note that in this expression we 



neglect the effects of nonlocal electrodynamics |27 which 



can not perturb the wave function g essentially in the 
limit e < KVa /S. when the minimum characteristic wave- 
length Imin ^ HVa/s exceeds the coherence length ^. One 
can see that our equations are analogous to the ones de- 
scribing the quantum mechanical motion of a massless 
Dirac particle with a charge |e| in the 'vector potential' 
8- = ~0oVx/(27r) of AB solenoids and the scalar poten- 
tial (p of 2D 'electrical dipoles' screened at a length scale 
Al- Both the solenoids and the dipoles are positioned at 
Vi. Each solenoid carries the flux quantum (j)Q. The ex- 
pression for a dipole moment reads: P = — 0.5?iyFyo/|e|. 
The Hamiltonian (|]) provides a simple tool for the study 
of 2D quantum mechanical effects in the QP motion and 
has been recently used as a starting point for both the 
analytical and numerical analysis of the band spectrum 
in regular vortex arrays with a rather small intervortex 
distance ii„ < Al Hj^H- 

Isotropic Dirac cone. Let us focus on the case of a 
single isolated vortex line. We start our analysis from 
the most simple isotropic limit Vp = Va = ^ (isotropic 
Dirac cone) and show that the low energy QP states near 
each node are strongly influenced by the pure gauge po- 
tential a due to the AB scenario. Indeed, if we neglect 



the potential of a 2D 'electric dipole' in (|i|), the scatter- 
ing cross section of a Dirac fermion in the AB potential 

appears to diverge for e ^ 0: — oc e^^ pSl-^. Intro- 

dO /■ — ' ' — " 
ducing a polar coordinate system (r, 9) with the origin at 

the vortex center and taking the order parameter phase 

in the form x = ^ one can obtain the eigenfunctions of 

the Hamiltonian Hq: 

oc laye^'^g:^ , (4) 

where Ji, is the ly-th Bessel function, m is an integer, 
\e\ = hVk, and L is the system size. QP states conflned 
near each node provide the following contribution to the 
local DOS: 



N = 



J2^J dk{\u^\^S{s-ek) + \um\'Sie + ek)) , (5) 



where Sk = iTiVk. The wavefunctions gln^ and gln^ with 
TO > are regular at the origin. The local DOS corre- 
sponding to a set of these regular solutions vanishes in the 
region r < fiV/\e\ and approaches the value Noo oc |e| for 
r ^ S-V^/kl {Noo is the DOS in the absence of vortices). 
The solutions with negative m diverge at the origin and 
are responsible for the formation of the nonzero DOS 
near the vortex inside the domain r < ?iy/|e|. Thus, the 
residual DOS Nq (at e = 0) is also determined by the con- 
tributions from the states with negative to. The crucial 
role in the formation of the residual DOS is played by the 
states ^f^i and g^l which appear to decay most slowly 
from the vortex center (as r^^/^) in the limit e ^ 0. Us- 
ing (^ one can obtain the following contribution to the 
residual local DOS: Nq ~ {hVr)~^. Let us emphasize 
that this contribution appears to be nonzero even in the 
large-r domain (r > Al) in sharp contrast to the behav- 
ior expected on the basis of the semiclassical model which 
takes account of the Doppler term \e\(p and neglects the 
AB effect. 

We now proceed with the analysis of the effect of the 
Doppler term in the Hamiltonian (^ on the behavior 
of wavefunctions in the small-r domain (r < Xp). For 
low energies e ^ hV/Xp the regular solutions (||,||) with 
TO > are only weakly influenced by the potential of the 
screened 'electric dipole'. On the contrary, for the wave 
functions with negative m (which diverge at r — > 0) 
the scalar potential ip can not be considered as a small 
perturbation and essentially modifies the solution in the 
domain r < Xp. In the low energy limit e <C TiV/Xl the 
solutions in this r-domain can be written in the form: 

5« = (1 + ^CT,)r^+l/2exp(-^0(l + a,)/2)G^(0) , (6) 

(7) 



3 



where the equation for G^{6) reads: 

d '- siTiO 
t-Q^Gf, + (1 + /i)(7^G^ + ^<7xG^ = . (8) 

A set of discrete quantum numbers fi is determined by 
the condition: G^{0) — G^{9 + 2ti). In the large-|/z| Umit 
the /i values can be calculated using the quasiclassical 
quantization rule: 



271- 



(1 + + ^siin?ed0 = 2?™ , 



(9) 



where n is an integer. The residual DOS near the vortex 
center is dominated by the states (^,^ with /i < —1/2, 
which are characterized by the power divergence at r — > 
0. The divergence should be cut off due to the matching 
with the solution inside the core, which results in a strong 
mixing of QP wavefunctions characterized by different /i 
values and corresponding to all four nodes. Outside the 
core the most slowly decaying wavefunctions correspond 
to the value /i = — 1 and can be found exactly using 
Eq.(|): 



-iO 



„(2) -ie-* 



(10) 
(11) 



where 7 = cos6'/2, and the constant C cx L~^l'^ is deter- 
mined from the matching with the solution in the large-r 
domain (r > A^). One can see that for intermediate dis- 
tances ^ <C r <C Ai these states provide a dominating 
contribution to the residual local DOS iVg ^ [fiVr)"^ , 
which coincides with the one predicted on the basis of 
the semiclassical model [|l^. Obviously the latter con- 
clusion is no more valid near the vortex core, where the 
admixture of the solutions decaying more rapidly from 
the vortex axis becomes substantial and can result in a 
narrow DOS peak similar to the one observed in high-Tj, 
compounds ]9|-pH. 

Despite of a possible formation of such a peak struc- 
ture the resulting low-energy states are not truly local- 
ized. The wavefunctions g'^-y and can be considered 
as leading terms in a large distance (r ^ ^) asymptotic 
expansion for the low-lying states and, thus, are respon- 
sible for the escape of QPs from the core. These conclu- 
sions can be supported by the analysis of the generalized 
inverse participation ratio, defined as p9|: 



/3cx 



(/(I^P 



For slowly decaying solutions glX and g'^i the /3 value 
appears to grow logarithmically as a function of an in- 
creasing system size L. Such logarithmic behavior is in 



a sharp contrast to the /3 oc divergence expected for 
localized states. The absence of the /3 cx scaling is in a 
good agreement with the numerical analysis |jl^ carried 
out in the limit L ^ A^, . 

Anisotropic Dirac cone. The physical picture sug- 
gested above can be generalized for the case of an 
anisotropic spectrum Vp 7^ Va- Such a generalization 
is of particular importance for the understanding of the 
electronic structure of vortices in high-Tc cuprates, where 
an estimate based on the results of thermal conductivity 



measurements |32 gives us a rather large anisotropy ra- 
tio: a = Vf/Va '^-^ 14. Taking an appropriate gauge 
transformation g = e^^ f one can choose the 'vector po- 
tential' a in the Hamiltonian ifg in the form: 



a = — 



00" (-yxo + xyo) 



2ti 



{ayf 



Using the scale transformation {ay — y, x ^ x) the 
Hamiltonian Hq can be reduced to the isotropic one 
with V ~ Vp. As a result, in the new coordinates 
the solution / in the large-r domain (where the Doppler 
shift is negligible) can be written in the form (HJ^) with 
k = |£|/(?iyF). The small distance asymptotic expan- 
sion can be analysed analogously to the isotropic case. 
The most slowly decaying solutions can be obtained from 
( |lo| , pl]) if we replace and r by 6' = tan^^{y/x) and 
r — y^x^ + y^, respectively, and take 



7 = 



2V^2— T 



tan ^(V"^ — IcosB) . 



Outside the core the expression for the residual local DOS 
(taking account of the contributions from all four nodes) 
reads: 



No 



+ 



V^x^ 



(12) 



The anisotropy of the Dirac cone results in the angular 
dependence of the local DOS which exhibits a fourfold 
symmetry with the maxima along the nodal directions 
(see Fig.l). 

Multi-quanta flux structures. The generalization of the 
above analysis to the case of a multi-quanta vortex with a 
winding number is straightforward. For a vortex car- 
rying an odd number of the flux quanta the effect of 
the AB potential on the local DOS outside the flux tube 
(r > Al) is the same as for a singly quantized vortex. 
On the contrary, a vortex with an even number does 
not cause the AB interference and, thus, can not pro- 
vide the slowly decaying contribution (iVp oc r~^) to the 
residual DOS discussed above. The most direct way to 
observe this odd-even effect is to consider a hollow cylin- 
der with a trapped magnetic flux $. The AB mechanism 
is the cause of an oscillating dependence Nq{^), which 
should be observable by any experimental technique sen- 
sitive to the residual DOS. The large odd-even effect in 
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the residual DOS is specific for superconductors with gap 
nodes and can be considered as a test for d-wave pair- 
ing. Obviously, the magnitude of such odd-even effect 
is strongly influenced by finite lifetime and temperature 
effects. One can expect that these mechanisms should 
suppress the slowly decaying contribution ( p^ ) to the lo- 
cal DOS in the domain r > TiVp / 'max\r ^ T], where F is a 
scattering rate. 

To sum up, we have discussed an extremely impor- 
tant role of pure gauge potentials in the behavior of low 
energy QP excitations in the mixed state of superconduc- 
tors with gap nodes. Within the BdG theory linearized 
in gradient terms we have considered the peculiarities of 
extended QP states for isolated flux lines. Such a model 
allowed to develop a 2D quantum mechanical description 
of these states, taking account of both the Doppler shift 
of the QP energy and the Aharonov-Bohm effect. It is 
hoped that the physical picture considered in this Letter 
can provide a starting point for the analysis of static and 
dynamic properties of the mixed state in various d-wave 
systems, including, probably, high- Tc copper oxides. 
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FIG. 1. The angular dependence of the residual DOS 
No{9)/No{0 = 0) in polar coordinates for a — 5 (dashed line) 
and Q = 20 (solid line). 
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